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When we consider a binary mixture of two chemical species 1 and 2, in equilibrium 
with its vapor, one of the problems, experimentally as well as theoretically, is to predict 
how the corresponding surface tension depends on the composition of the mixture. Some 
relationship is expected which would give this surface tension, as an interpolation between 
the surface tensions of the two species when they are chemically pure. In this note, we 
briefly discuss various relation-ships based on thermodynamical considerations as well as 
other ones obtained more recently in the frame of Solid-On-Solid models of interfaces 
and bulk statistical mechanical models of binary mixtures [2] ■ 

Thermodynamical semi-empirical equations 

Let us use T(i,2)|o to denote the surface tension of a mixture of two species 1 and 2 
and let Ti| and r 2 |o be the surface tensions of each species. In this section, we present 
several equations which have been derived, according to different assumptions, by using 
thermodynamical considerations. 

For ideal or nearly ideal solutions, a fairly simple treatment, due to Guggenheim [2], 
leads to the following equation 

e -0aT (1 , 2) | O = Cie -paT M0 + C2e -/3ar 2 | ^ 

where c\ is the fixed molar fraction of species 1 in the (1,2) mixture, c 2 = 1 — c%, the 
fixed molar fraction of species 2, a is the mean surface area per molecule, and (5 = 1/kT 
is the inverse temperature. 

A very simple relationship for the so called regular solutions comes from Prigogine and 
Defay [3], who proposed the equation 

T (l,2)|0 = c l T l\0 + C 2T 2 \0 ~ KciC 2 (2) 

with K a semiempirical constant. 
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A simple treatment due to Eberhart |2j assumes that the surface tension of a binary 
solution is linear in the surface composition, that is 

T (l,2)|0 = c l r l|0 + C 2 T 2|0 ( 3 ) 

where the cf , % = 1,2, denote the mole fraction near the surface of phase separation, and 
that the ratio c\/c\ is proportional to the ratio c|/c2- 

Finally, when the surface tensions t\\q and T2\o differ appreciably, a semiempirical 
equation attributed to Szyszkowsky (j^j, jjj) gives: 

T(1 < 2)|0 -1-I?ln(l + ^) (4) 



7i| \ A. 

where two characteristic constants A and B of the compounds have been used, and C2 is 
the concentration of the species with the smaller surface tension. 

We refer the reader to Adamson's book |B] (Chapter III, Section 4), and references 
therein, for a detailed discussion of the above equations. 

Let us also mention that an extensive development for various types of non ideal 
solutions has been made by Defay, Prigogine and co-workers: see the monography 0|. 

Solid-On-Solid (SOS) models 

Let us first consider the SOS model in dimension d = 2. We let ho, ■ ■ ■ , /jjv £ K be a 
collection of heights describing an interface. For simplicity, we assume that the energetic 
cost of the interface is proportional to its length, namely given by 

N-l 

H(J a ; h --- h N ) = J a 22 t 1 + \ h i+l ~ hi\] (5) 
i=0 

where J a represents the energetic cost per unit length for the interface. The associated 
density of free-energy or interfacial tension is then defined by 

r(PJa) = ~\ Km -j- In Z(N, (3J a ) (6) 

p N—>oo iv 

where Z{N,(3J a ) is the following partition function 

/+oo r+oc 
dho--- / dh N e-P H{Ja > h °- hN) 5{ho - 0)5(h N - 0) 
-oo J — oo 

The quantities r(/3Ji) and TifiJ?) will represent the surface tensions of two species 1, and 
2. 

To model the interface for binary mixtures, we will use the disordered generalization 
of the SOS model. Namely, we consider that the coupling per unit length is no more a 
constant but a random variable J. This variable may take two values J\ and J2 with 
probabilities c\ and C2 = 1 — c\ , where c\ is physically the mole fraction of particles 1 in 
the bulk of the (1,2) mixture. This allows us to represent the interface by two sets of 
independent random variables {ho, . . . , h^} and { J°, . . . , J N_1 }. In this way, the interface 
can adjust its height h and, by moving or not molecules, also the corresponding energetic 
cost J taking into account the fixed concentration. In this approach the molecules of 
the mixture are not distinguishable and have similar size. Each site i is occupied by one 
molecule. The energetic cost of this interface is given by the Hamiltonian 

N-l 

H(J° ■ ■ ■ J N ~ 1 ; h ---h N ) = J2 J*[l + \hi+l ~ tH\] (7) 

i=0 
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where J l may take the value J\ or J 2 with probability c\ or Cq. = 1 — c±. The associated 
partition function is given by 



+00 r+00 



00 



Z N (J° • • • J"' 1 ) = dho--- dh N e-P H ( j0 --- jN - 1 -' h °--- h ^8(h )6(h N ) 



It is well known that random systems are often related to disordered systems for which 
one has introduced the notion of quenched and annealed disorder. For the annealed case, 
the couplings are considered to be random and will be treated in the same way than the 
heights. For the quenched disorder, the couplings are frozen in a given configuration. 
There are then two ways to define the associated free energy density 



and 



-™o led = ^ MZn(J° ■ ■ ■ J N ~ 1 )) (9) 



where the average ( • ) has to be taken with respect to the coupling distribution. 

Using quenched disorder to compute the free energy, we obtain the following equation 
to express the surface tension of the mixture like a convex combination of the pure 
component surface tensions 

quenched = ^(fij^ + c 2 t((3J 2 ) 

This kind of formula is therefore valid when the system finds its equilibrium position 
within the configurations for a given set of couplings. 

The other approach, i.e. if we use the annealed disorder, leads to : 

e -^— = cie -^. h) + C2e - Mf sj 2) 

Here, the variables ,P and hi are treated on an equal basis. This implies that the molecules 
at the interface are sufficiently mobile to allow the interface to adjust itself, both in heights 
and in composition. 

Bulk model of binary mixture 

This Section is devoted to a discussion of the problem within a lattice bulk statistical 
mechanical model describing the binary mixture in equilibrium with its vapor. Studies 
of various models of binary lattice gases can be found in Refs. jlUl lllj . 

Here, we consider a lattice gas system with two kinds of particles, where each lattice 
site can be in one of the three states, 0, 1,2, interpreted, respectively, as an empty site, a 
site occupied by a particle of the first kind of the model, and a site occupied by a particle 
of the second kind. Whenever the particles 2 are not allowed the system reduces to the 
usual Ising model, in its lattice gas version, with coupling constant Ji/2, and analogously, 
when particles 1 are not allowed, it reduces to the Ising model with coupling constant 
J2/2. Namely, to each site x € d = 2,3, we associate a variable s x taking values in 
the set = {0, 1,2}. We will say that the site x is empty when s x = and that it is 
occupied otherwise. The energy of a configuration sa = {s x } x& \ in a finite box A C TL d 
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is defined by 

H\(s A ) = E(s x ,s y ) 

{x,y)cA 

E(s x , s y ) = Ji [8(s x , l)5(s y ,0) + 0)<5(s y , 1)] 
+ J 2 [5( aa ., 2)<5(s y ,0) + 5{s x ,0)5{sy, 2)] 

where denote nearest neighbour pairs, 5 is the usual Kronecker symbol S(s,s') = 1 

if s = s' and 5(s,s') = otherwise; Ji and J2 are positive constants. Notice that the 
bond energy satisfies: 

E(0, 0) = £7(1, 1) = £7(2, 2) = E(l, 2) = E(2, 1) = 
£7(0, 1) = £7(1,0) = Ji,£7(0,2) = £7(2,0) = J 2 

Fixed densities of the three species are introduced through the canonical Gibbs ensem- 
ble of configurations sa such that 

5^(^,0) = No, '£ / S(s x ,l) = N 1 , J2 S ( S *> 2 ) = N 2 

x&A xeA xgA 

where the sum iVo + N\ + N 2 equals the number of sites of A. The associated partition 
functions with boundary condition be are given by 



£ e-Ws^^ (si)1)|]Vl U^ ( 
i A en A VxgA / \xeA 



ZUA;N U N 2 ) = e -^A(s Aj(5 1 ) ^ (Sx)1))iVl ) 5 [ > ^( S:E ,2),iV 2 ) x &c (s A ) 

where x 6c (sa) is a characteristic function standing for the boundary condition be. We 
shall be interested in particular to the following boundary conditions: 

• the empty boundary condition: x° mp ( s A) = FLe^A 8{s x , 0) 

• the mixture boundary condition: x mixt ( s A) = ELeaA-U ~~ 8(s x ,0)) 

• the free boundary condition: x & ( s a) = 1 

where hereafter, the boundary dA of the box A is the set of sites of A that have a nearest 
neighbour in A c = Z d \ A. 

The free energy per site corresponding to the above ensemble as a function of the 
densities p\ and p 2 of the particles 1 and 2 is 

f( Pl ,p 2 ) = lim-^-lnZ bc (^[pi\Ml[P2\M\) (10) 
A|Z d p|A| 

where [ • ] denotes the integer part and the thermodynamic limit A f Z rf is taken in the 
van Hoove sense [12] . 

It is convenient to consider also a grand canonical Gibbs ensemble, which is conjugate 
to the previous ensemble, and whose partition function, in the box A is given by 

E bc {A;p llf i 2 ) = Yl e-^^'^^^Ai^.iJ+wE^A 5 (^> 2 ) (11) 

s A GSl A 

where the real numbers p\ and p 2 replace as thermodynamic parameters the densities p\ 
and p 2 . The corresponding specific free energy, the pressure, is the limit 

p(Mi>M2)= lim j T - T lnH 6c (A;// 1 ,/x 2 ) (12) 

A\7L d |A| 
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Limits and (|12j) . which define the above free energies, exist. They are convex func- 
tions of their parameters and are related by the Legendre transformations 

p(Vi, M2) = sup [mpi + P2P2 - /3f(pi,P2)] (13) 

0f(pi,P2) = SUp [pipi + p 2 P2 ~ p{pi,P2)} (14) 

Finally we introduce the finite volume Gibbs measures (a specification) associated with 
the second ensemble: 

e -/3H A (s A ) bc( \ 

< c (sa) = p X SA) (15) 

where H\(s A ) = H A (s A ) - ^ Y,xeA $( s x> l ) ~ *f XLgA <H s z, 2 )- Tne y determine by the 
Dobrushin-Landford-Ruelle equations the set of Gibbs states Qp{7i) on 7L d corresponding 
to the Hamiltonian H at inverse temperature (3 (see e.g. 1131). A Gibbs state P £ Opi^i) 
which equal the limit lim A | Z d P A (-), is called Gibbs state with boundary condition b.c. 

Ground states and low temperature analysis 

In the zero temperature limit the Gibbs state with empty boundary condition is con- 
centrated on the configuration with empty sites: 

lim P A c (emp A ) = 1 (16) 

where emp A is the configuration where all the sites of A are empty, and this limit vanishes 
for any other configuration. Gibbs states at f3 = 00 are called ground states. 

Let, i? A ixt (c) = {s£ f2 : \/x £ A, s x /=0}, be the restricted ensemble of configurations 
in A with non empty sites, and i? A lxt (c), < c < 1 the subset of configurations of i? A lxt 
with exactly [c|A|] = N sites occupied by a particle of the specie 1 (and |A| — [c|A| sites 
occupied by a particle of the specie 2). 

With the mixture boundary conditions one gets by by Stirling's approximation and 
the principle of maximal term that 

lim Wt*{Rf*{c)) — 1 (17) 

for c = gpf^Ma ■ This means that the ground state with mixt boundary conditions is 
concentrated on the restricted ensemble i? mixt (c) of configurations of non empty sites 
with concentration c of particles 1 and concentration 1 — c of particles 2. 
With free boundary conditions, one has for c = eM1 e ^ 1 gM2 

lim p{ r (emp A ) = 1/2, lim p{ r (Rf xt (c)) — ► 1/2 

Thus, with the above considerations, we get that for e Ml +e^ 2 = 1, the configuration with 
empty sites coexists with the restricted ensemble itl mixt (c). 

This analysis can be extend to the Gibbs states at low temperatures by using Pirogov- 
Sinai theory [Tl|. Actually, this theory allows to show that the low temperature phase 
diagram of the model is a small perturbation of the diagram of ground states the coexis- 
tence line given by the equation 

l n ( e Mi + e *4) = e /4-2<Wi + 6 a4-2<W2 



where J = min{ J\, J2}. 
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In particular, introducing the infinite volume expectation ( • ) (fii, H2) associated to 
the Gibbs measure (fTKj) . ( • ) bc (/xi, ^2) = nm ATZ d J2s A en A ' ^a c ( s a)i we have for any t > 0: 

(5(s x , 1) + 5( fla! , 2)> mixt ( / x 1 + 1, ,4 + 1) > 1 - O (e~ 2d ^ 

(«5( S:c ,l) + 5( S:E ,2)) em P(^-t,^-t) < o(e" 2 ^ J ) 

showing that the model exhibits at low temperature a first order phase transition at the 
coexistence line where the pressure is discontinuous. 

Surface tensions 

To introduce the surface tension between the mixture and the vapor, we consider the 
parallelepipedic box: 

V = V L , M = {(zi, x d ) e Z d : |xi| < L, i = 1, d - 1; — M < x d < M - l} 

and let d + V (respectively d-V) be the set of sites of dV with Xd > (respectively < 0). 
The boundary condition, x mixt > cm P(s v ) = Y[ xedV (l - S(s x , 0)) ]J xea+v $( s xi °)> enforces 
the existence of an interface between the mixture and the vapor and the interfacial tension 
between the mixture and the vapor is defined by the limit 

1 ,. 1 , ^mixt,emp(^; M*, ^2) 



r(l,2)|0 = -- lim n — ^r~r In 



/3i^oc,M^oc (2L + l) rf -! (Z^(V;u$,u*)Z emp (V;ut,u* 2 )y/ 2 

As mentioned previously, whenever either the particles 1 or the particles 2 are not allowed 
the system reduces to the usual Ising model in its lattice gas version. Thus, to define 
the surface tensions between each species of the mixture and the vapor, we introduce the 
configurations ny G {0, 1}^ of the lattice gas and the following partition functions 

Q a (V) = ^ e^^^c^ 1 -™^^ 1 -™^ 1 Yl n x 
n v e{o,i} v x &dv 

n A e{0,l} v x&d-V xed+V 

for a = 1 and a = 2. The interfacial tension between the species a = 1,2, and the vapor 
is the limit (H5HE|) 

T n n = hm lim . , z 

°' U L^ooM^oo (2L + l)^ 1 

where F a (V) = — ^ In ^4yy ■ It is well known that the ratio Q a ,o{V)/Q a {V) can be 
expressed as a sum over interfaces which in this case are connected set of bonds or 
plaquettes of the dual lattice |17| I18j . Extracting the energy of the flat interface, the 
system can be written as a gas of excitations leading to F a (V) = J a (2L + l) d ^ 1 + F£ K (V). 
In two dimensions, is the free energy of the gas of jumps of the Gallavotti's line |17j . 
In three dimensions, is the free energy of the gas of walls of the Dobrushin's interface 
|18j . In both cases these free energies can be analyzed by cluster expansion techniques at 
low temperatures. Namely, the specific free energies T a = lim^oo F^ xt (V)/(2L + l) d_1 
exist and are given by convergent expansions in term of the activities e~^ Ja , giving 

T a ,0 = Ja + Fa. (18) 
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In addition 

- PF a = 2e~P Ja + 0(e~ 2l3Ja ) for d = 2 (19) 
-#F Q = 2e- 4/3Ja +0(e- 6/3Ja ) for d = 3 (20) 

Furthermore, in two dimensions the surface tension defined above is known to coincide 
with the one computed by Onsager We thus have an exact expression for t Q) o, and 
for J- a , namely, f3J- a = lntanh(/3J a /2)for (3J a larger than the critical value ln(l + \/2). 

Main result 

The above surface tensions Tn,2)|o> r i,o and T2,o, are proved to satisfy, whenever (3 is 
large enough, the equation 

e -0(T(i,a)|o--n = c * e -^(^|o-^i) + c^-^^k,-^) (21) 

The quantity T is the specific free energy (which can be expressed as a convergent series at 
low temperatures) of a gas of some geometrical objects called aggregates. In dimension 
d = 2, those aggregates are the natural generalizations to our model of the jumps of 
Gallavotti's line and the leading term of the series giving this free energy T is 

2 C * e -Wi + c * 2e -WJ2 

~P d[e-P J i +c* 2 e~P J 2 

In dimension d = 3, they are the natural generalizations of the walls of the Dobrushin's 
interface and then the leading term of the series is 

1 c\e-^ h + c* 2 e-^- h 1 {c\e~ 2 ^ + c\e~ 2 ^f 

The coefficients c* and c 2 are related to the concentrations c\ and c 2 of the particles 1 
and the particles 2 through the equation 

c * = e f?-PO*I.f5) } i = 1,2 
This equation gives at low temperatures: 

1 - ( Cl e-^ + c 2e -^) M - JWcie-^ ( Cl e^ Jl + c^ 2 ) 2 "" 1 
-2(d + l)a (de-^ 1 + C2 e-^) 2d + O ( e -Vd+W™*{Ji,J^] (22 ) 

References 

[1] J. De Coninck and J. Ruiz. Interfacial tensions for binary mixtures versus Solid-On-Solid models. 
Preprint, 2003. 

[2] J. De Coninck, S. Miracle-Sole, and J. Ruiz. On the statistical mechanics and surface tensions 
of binary mixtures. Preprint CPT-Marseille-2004/P-015, mparc 04-206, submitted to Journal of 
Statistical Physics. 

[3] E. A. Guggenheim. Trans. Faraday Soc, 41:150, 1945. 

[4] R. Defay and I. Prigogine. Trans. Faraday Soc, 46:199, 1950. 

[5] J. G. Eberhart. J. Phys. Chem., 70:1183, 1966. 

[6] B. von Szyszkowsky. J. Phys. Chem., 64:385, 1908. 

[7] H. P. Meissner and A. S. Michaels. Ind. Eng. Chem., 42:2782, 1949. 

[8] A. W. Adamson. Physical Chemistry of Surfaces. John Wiley and Sons, New- York, 1997. 

[9] R. Defay, I. Prigogine, A. Bellmans, and D. H Everett. Surface tensions and absorption. Longmans, 

Green and Co., London, 1966. 
[10] J. C. Wheeler and B. Widom. J. Chem. Phys., 52:5334, 1970. 



8 



JEAN RUIZ 



[11] J. L. Lebowitz and G. Gallavotti. Phases transitions in binary lattice gases. J. Math. Phys., 12:1129, 
1971. 

[12] D. Ruelle. Statistical Mechanics: Rigorous Results. Benjamin, New- York, Amsterdam, 1969. 
[13] D. Ruelle. Thermodynamic Formalism. Addison Wesley, Reading, 1978. 

[14] Ya. G. Sinai. Theory of Phase Transitions: Rigorous Results. Pergamon Press, London, 1982. 

[15] A. Martin LI, G. Gallavotti, and S. Miracle-Sole. Some problem connected with the coexistence 

of phase in the Ising model. In Statistical mechanics and mathematical problems, Lecture Notes in 

Physics, volume 20, page 162. Springer, Berlin, 1973. 
[16] J. Bricmont, J. L. Lebowitz, and C.-E. Pfister. On the surface tensions of lattice systems. Ann. N. 

Y. Acad. Set., 337:214, 1980. 
[17] G. Gallavotti. Phase separation line in the two-dimensional Ising model. Commun. Math. Phys., 

27:103, 1972. 

[18] R. L. Dobrushin. Gibbs states describing the coexistence of phases for a three-dimensional Ising 
model. Theory Prob. Appl, 17:582, 1972. 



